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We report on a construction for smectic blue phases, which have quasi-long range smectic trans- 
lational order as well as three dimensional crystalline order. Our proposed structures fill space by 
adding layers on top of a minimal surface, introducing either curvature or edge defects as necessary. 
. . . We find that for the right range of material parameters, the favorable saddle-splay energy of these 

CO I structures can stabilize them against uniform layered structures. We also consider the nature of 

, curvature frustration between mean curvature and saddle-splay. 

o 
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I. INTRODUCTION 



Topological defects and liquid crystal phases are hopelessly intertwined. Historically, the nematic phase derived it's 
' name from the observation of long disclination lines. In some liquid crystal phases, such as the twist-grain-boundary 
O , (TGB) phases and blue phases, these defects are not simply artifact or nuisance, but instead act to stabilize the 
director field configuration. While the TGB phases have smectic order, the blue phases have purely nematic-like order 
' created by a long range, triply periodic lattice of line defects. Two of the nematic blue phases possess cubic symmetry 
. {BPl and BP2) while the third {BP3) is thought to be an isotropic melt of double- twist cylinders Recently, new 

^ ' phases of matter have been identified that possess the quasi-long range translational order of smectics in addition 
to three-dimensional orientational order. These phases, dubbed the "smectic blue phases," have been observed for 
^ molecules in the chiral series FH/FH/HH-nBTMHC, where n is the aliphatic chain length. Three distinct smectic blue 
O phases have been observed near the isotropic transition of these compounds: BPsmA^ has cubic symmetry, BPsm'2 
has orthorhombic symmetry and BPsm'i is isotropic 13. The precise physical properties of these materials have been 



the study of intense investigation in recent years 0, Q| • 
T-H , In general, since smectic order is incompatible with cubic symmetry, it is expected that any triply periodic crystalline 
^ ■ structure must include smectic dislocations as well as disclinations. However, attempts to construct smectic double- 
2^ . twist cylinders @ and assemble them into traditional blue phase structures present a variety of difficulties, most 
^ ^ ' notably a disagreement with precise experimental details ^^l. In previous work we proposed a new model for the 

. smectic blue phases. Our construction filled space with concentric minimal surfaces wrapping a lattice of intersecting 

line defects. We found that when the saddle-splay constant was large enough, these new structures were stable. 
Though the new materials were chiral, our construction did not rely on macroscopic chirality as in the traditional 
: blue phases. Instead, the smectic compression and bending energies set the length scale of our solutions. 

Here, we refine our earlier model through variations on the original construction. We find that by allowing edge 
dislocations our phase is stabilized for even small (negative) values of the saddle-splay constant K24- This article 
is organized as follows. In Section ^] we discuss the rotationally-invariant energetics of layered systems and derive 
detailed equations for the geometric frustration between curvature and uniform layer spacing. Next, in Section UTTI we 
calculate the energy and stability of likely smectic structures based on our original construction [lOj . To supplement 
Q our analytical calculations, we present the results of simplified numerical solutions for the three dimensional smectic 
O structure, which explore relaxation of the smectic layers away from our constructions. In Section llVl we present a new 
^ ■ construction for filling space based on uniform layer spacing away from a minimal surface. First we derive formulae 
• I— I for the evolution of curvature fields in layered space and then we employ the Weierstrafi analytic representation of 
minimal surfaces to calculate the energies for our proposed phase. We complete our description by calculating the 
^ , core energy and argue that tilt-grain boundaries form at the cores of our structures. Numerics are also presented for 
■ - - ■ this construction. Finally, in Section^we calculate the Fourier transform of the smectic density and compare it with 
experimental X-ray results. In Appendix A we derive curvature evolution equations in curved space and in Appendix 
B we review the Weirstraf3 representation for completeness. 



II. ENERGETICS AND CONSTRUCTION 



The similarity between the crystal structures and phase diagrams of the nematic and smectic blue phases suggests 
that, like the nematic blue phases, the smectic phases are stabilized through saddle-splay. The key to our construction 
is the observation that saddle-splay and Gaussian curvature are identical 11 1| for layered systems with uniform spacing. 
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The saddle-splay energy of a unit director field N is |3] 

i^ss = ^24 /d^xV-KN- V)N-N(V-N)], (1) 



where if 24 is a Frank constant. In a layered system, we can rewrite the expression for _Fss in ^ more useful form 
by employing a local coordinate system where one direction is parallel to the local layer normal, N. This frame is 
appropriate in the limit where the nematic director and the layer normal are locked. Then Eq. ^ becomes 



Fss = -'2K24: dn dxdy^/gj^x~y) an{x,y)Kn{x,y) (2) 



where n is the Lagrangian coordinate which labels the layers, an{x,y) is the local layer spacing at {x,y), Kn is the 
Gaussian curvature of the n^^ surface, and gn is the determinant of the two-dimensional, induced, surface metric. In 
the special case that a„(a;, y) is constant, the integral becomes purely topological: a consequence of the Gauss-Bonnet 
theorem is that for a surface of genus g the integrated Gaussian curvature is 47r(l — g). Since 5 > for any 
infinite surface [l2j |. i^gg is large and negative when K24 < and the unit cells contain surfaces with large genus, 
i.e. many handles and holes. Note that here, the saddle-splay is a measure of the layer normals and not the nematic 
director. When the nematic director follows the layer normal these are, of course, equivalent. As discussed in |0, in 
type-II smectics it is possible for the saddle-splay of the director field to differ in its precise numerical value from the 
saddle-splay in the layers. 

In addition to the saddle-splay, we must include the rotationally invariant bulk free energy: 



Fsm ^ I d^'x \j 



+2KiH'\ (3) 
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where the smectic density is p esc cos (27r$/ao), ^{x, y, z) is a phase field, ag is the layer spacing, B is the compression 
modulus, Ki is the bend modulus, and H — ^V-N is the mean curvature of the layers. If the two principle curvatures 
are Ki and K2, Kg — K1K2 and H ~ [ki + K2)/2 [TiI |. 

Together, Eqs. 10) and Q favor configurations with uniform layer spacing, low mean curvature and high Gaussian 
curvature. However, there is an unavoidable geometric frustration between these three terms. When a surface is 
displaced along its normal by (5a, the changes in the metric and curvature tensors are |13|: 

5gij — 2KijSa 

Snij — —ViV jSa + Kikn'^jSa (4) 
It is straightforward to find the variation in the mean and Gaussian curvatures: 

5H = [KG~2H^]5a-^g'^V^VJSa 
SKg = -2KGHSa~ k'^V.VjSa. (5) 

where Vi is the covariant derivative on the surface and k*^ = e^'^e^'Kki/\g\, where e*-' is the completely antisymmetric 
tensor. 

If we impose smectic order, then a may also be interpreted as the normal distance between adjacent layers. When 
the spacing is uniform, N • V$ is constant and Wi6a = 0. In this case, the first result in Eq. (0) implies dH/da = 
if and only if the principal curvatures Hi satisfy 

kI + kI=0 (6) 

i.e. Ki = K2 = 0. Thus if Kg =/= and the layer spacing is uniform, then dH/da =/= 0; Gaussian curvature leads to 
mean curvature. 

For uniform spacing the evolution equations become particularly simple. The variation Sa is constant and the 
equations can be integrated: 

H + aKG 

Hia] — 

^ ' l + 2aH + a^KG 



FIG. 1: Minimal surface repeat units: (a) P surface, (b) I-Wp surface. 



Alternatively, this evolution follows from the observation that if the principle radii of curvature are i?i = and 
i?2 = for one surface then the radii for the surface displaced by a along the local normal are Ri(a) = Ri + a. 
Thus, in layered systems, Gaussian curvature in one region implies mean curvature in another. Furthermore, as we 
continue to develop the initial layer, there will be a curvature singularity at a distance a = {—H ± V H'^ — Kg) /Kg 
normal to the original surface, i.e. at one of the radii of curvature. 

With these geometric constraints in mind, we construct a solution which strikes a promising balance between mean 
and Gaussian curvatures. To do this, we build the smectic order by adding layers on top of a triply-periodic minimal 
surface. By starting with a minimal surface, we bias our structures to have a low total mean curvature energy. In 
we based our construction on the Schwartz P surface, pictured in Fig. ^a). In this paper, we present calculations 
for both the P surface and Schoen's I-Wp surface, pictured in Fig. ^b), a surface with a larger genus and thus a 
lower saddle-splay energy. The P surface is an archetypal triply-periodic minimal surface and will provide a basis of 
comparison with earlier results, while the I-Wp surface shares the symmetries of experimental Smsp systems the 
latter will also prove to be more stable than the P surface. We will consider two variations of this construction: in the 
first we will allow the layer spacing to vary, while in the second we do not. In the former we find that the cores are 
composed of topological line defects, while in the latter we find that domain walls form in the relatively larger cores. 

III. SMALL CORE MODEL: SELF SIMILAR LAYERING 

A. Geometric Construction 

In this section we consider space-filling, layered structures in which each layer is continuous and every layer has 
the same global topology. We begin with a minimal surface and devise an explicit construction for filling the region 
away from that surface. We note that our bulk phase must contain line defects, since at some point the curvature 
will diverge, as is implied by Eqs. Q). By design, the layers immediately surrounding these line defects will have the 
same topology as the original surface. Therefore, the next logical step in our construction is to guess the optimal line 
defect structure. 

Besides selecting the correct topology, we should choose a defect complexion that minimizes the compression energy 
in the region between the minimal surface and the defects: it should, as much as possible, be equidistant from the 
minimal surface at each point. Thus, an ansatz for the optimal structure is generated by uniformly translating the 
minimal surface along its normal until it self-intersects, as shown in Fig. |21 The line defects will be disclinations of 
charge -1-1 or -1-1/2. In the case of charge -1-1 disclinations the layers immediately surrounding them are tight cylinders. 
From Fig. 121 a), we see that the defects in the P surface are all -1-1. From Fig.[2Ib) and (c), we see that the defects 
inside the I-Wp surface should have charge -1-1, but that outside, the translated layer approaches itself in a plane 
rather than on a line and so the disclination charge is -1-1/2. Fig.|2Id) and (e) show our proposed defect structures. 
For the defect structure outside the I-Wp surface, the location of the intersection of defect lines is a variable chosen to 
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(a) 




(d) 




FIG. 2: In (a) we show the P surface translated inward until self intersection. In (b) and (c) The I-Wp surface is translated 
inward and outward, respectively. Figures (d) and (e) indicate the line defect structure for our construction, based on skeletal 
graphs. 

minimize the total compression energy. It should be noted that in the case where all defects happen to be charge +1, 
our requirements of topology and equal spacing are exactly those best satisfied by the skeletal graph of the minimal 
surface 0|. We take the core region around the defect lines to have a radius of order ag, the smectic layer spacing. 
We model the core as a melted region of smectic with a free energy arising from the condensation energy. 

Finally, we fill the region between minimal surface and defects with continuous layers. We desire to have as little 
mean curvature as possible in this region, so we fill the region with dilated copies of the original minimal surface. 
In order to do this, we must cut the continuous minimal surface into smaller patches: we cut the surface along lines 
where one surface tangent is parallel to the line defect structure, and dilate it so that it shrinks onto the vertices where 
defect lines meet. The separated minimal patches are connected together by cylindrical patches parallel to the line 
defects. When there are disclinations, we must also add some flat patches to complete the surfaces. Thus the surface 
tangents are continuous everywhere, and the intermediary structures transform smoothly from the minimal surface 
on the outside to a set of cylinders surrounding line defects on the inside. This construction fills space completely, 
at the expense of uniform spacing. Moreover, we have filled a large part of space with minimal surfaces and so the 
curvature energy is likely to be as small as possible, given the necessary curvature defects. 

Fig. 13 shows how this construction works for the P and I-Wp surfaces. For the P surface it is particularly simple; 
its construction is show in Fig. |3Ia). The line defects are along the edges of the unit cell and through the center 
parallel to the edges. The inside of the P surface is filled with dilated copies of the P surface repeat unit, cut at its 
intersection with the walls of the unit cell. The edges of the dilated P surface copies are connected to neighboring 
cells with cylinders. The outside of the P surface is a translated copy of the inside. The I-Wp surface is more 
complicated, requiring different constructions for outside and inside as well as having both +1 and +1/2 defects. 
Fig. I2b) illustrates the construction inside the I-Wp surface: the minimal surface is divided at cuts that can be 
smoothly attached to cylinders, then dilated towards either the center or corners of the unit cell. Fig.|3Ic) shows the 
construction outside the I-Wp surface. The surface is cut along its mirror symmetry planes, x = y, x = — y, y = z 
etc. and x = 0, y = and z = 0. The resulting patches are dilated around the proper vertices in Fig.jSJe). The edges 
of the patches are connected together with cylinders parallel to the -1-1/2 edge defects. Finally, the free edges of the 
cylinder are connected together into planes. 

Most of the energetics of this structure are simple to calculate; exact values are tabulated in TableQ] It is instructive 
to derive approximate expressions for these energies so that we may quickly assess the feasibility of a given minimal 




FIG. 3: Small core construction: (a) Partial layers inside the P surface, (b) One layer inside the I-Wp surface, and (c) Partial 
layer around one octant outside the I-Wp surface. Each color corresponds to a different connected layer. In (b) and (c), the 
boundaries between minimal surface patches, cylinders and planes are drawn in heavy black. The defect networks and vertex 
points are also shown. 
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TABLE I: Smectic energy per unit cell for P and I-Wp surface constructions. The first column gives the topological genus per 
unit cell of the structure. Fb is the total curvature energy, Fss is the saddle-splay, Fcore is energy of the defect cores, and Fc 
is the total compression energy. 

surface as a starting point for our construction. 

The total curvature energy comes from the contribution of the cylindrical connecting patches. It is roughly 
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where the sum is over defect lines, li is the length of line i, Ci is the charge of the defect, pc is the core diameter of 
the line defect, and (R) is the average radius of cylindrical regions. The factor oi 1{R — r) / R above arises from the 
dependence of the cylinder length integration on the radius r. The core energy, vifhich is simply a condensation energy 
density, is approximately: 
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where e is a line tension. The first term is the energy of the defect lines, the second is the energy of the central, 
melted region around vertices in the defect network. 

Since the saddle-splay is a surface term, we calculate it by integrating the strength of the nematic defects along 
the total length of defect core. Because the saddle-splay is a total divergence, the volume integral becomes a surface 
integral around the defect: 



Fss = 2i^24 / HdS 
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(10) 



Comparison to Eq. shows that these smectic structures will favor charge +1/2 structures when possible. Also, for 
all constructions which only utilize -f 1 defects, the ratio of saddle-splay to curvature energy will be approximately 
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The only energy not explicitly determined by our construction is the compression energy. Though the layering 
is largely constrained, there is still one free variable which gives the relative dilation of consecutive surfaces. We 
minimize the compression energy with respect to the relative smectic layer spacing as a function of position. In the 
following we show how we can take advantage of the dilational symmetry of our construction to greatly simplify the 
exact calculation of optimal compression energy. 

First we consider the compression energy in a subregion which is filled with self-similar, radially dilated patches of 
minimal surface. Each successive layer in this region is a smaller version of the last, and all are centered on a common 
origin. Thus, if the radial coordinates of the outermost surface of this region are specified by a function ro(0, 0), then 
the radial coordinates of any interior surface are given by r = ('ro(0, 0), where C e [0,1] parameterizes the layers. 
Furthermore, every point within this region has a unique value of ^, </>, and ^, and conversely the values of 0, 0, and C 
uniquely specify the position of any point in the region, so we can use these three variables as a coordinate frame on 
the patch. In these coordinates, each layer is a surface of constant C. 

In the continuum description of the smectic, the layers are surfaces of constant $. Since $ is constant on each 
smectic layer, which is, in turn, a layer of constant ^, it follows that <I> = is purely a function of Our goal is to 
find the form of ^(C) which minimizes the compression energy. This is non-trivial, since the compression energy goes 
as V$, which is not purely a function of However, we can find a combination of |V3>| times a function of angle 
which together is constant over each smectic layer. First we observe that 

$(C) - <& 



-|*(C)-o(M)^$(;^) (11) 

Then we note that for level surfaces of $ the field of unit normals is N = V<I>/| V$| to write: 

(r • N) |V$| = r— $ ( 1 = — $(C) C— $(0 (12) 

This shows that the combination (r • N) | V^j depends only on (, and so is constant over each smectic layer. We can 
therefore write 



{V^f =p{9,4>)A{0, (13) 



with 



ro(go,(/.o).N(go,0o) ^' 
ro(0,0)-N(0,0) . 

l2 



(14) 



A(C) = An[V^Y , (15) 

where (^Oi'/'o) is a reference direction and An is a normalization constant which adjusts for differences between 
reference directions on different patches. The absolute value of An is unimportant, but its relative value on different 
surface patches must be chosen for a consistent definition of A(C) across the entire smectic surface. Thus we have 
separated the angular dependence out of the volume integral of V$. The result is 



I^C patch 2 " 



J dC (/o - 2/iA(C) + hA^C))] , (16) 
where the In are moments of the minimal surface shape, 

In^ f dnrl{0,cp)p^{e,c^) 



Simple consideration shows that a similar separation occurs in regions with cylindrical symmetry, in which case 
r = C^o(0) and |V$| is independent of the local z coordinate. The equations are similar to the above. Numerical 
values were calculated with the aid of the the Surface Evolver software package 15]. We minimize the total Fc by 
varying with respect to A(C). The results are tabulated in Table The total compression energy had only a very 
weak dependence on the core size ao for ao <C i. 

Calculations of the compression energy depend greatly on the particulars of the surface shape, so we cannot give an 
accurate approximate expression for Fc- In general, however, the total compression energy should increase rapidly 
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FIG. 4: Stability diagram for unit cells with L = 50ao. Our constructions based on P and I-Wp surfaces can be stabilized against 
flat layering in the regions of parameter space to the upper right of the respectively labeled lines. The quantity A = y'Ki/B 
is the penetration length. 



with genus. Since all layers in this construction have the same topology, there must be just as many layers in the 
cylindrical regions, of radius (R), as there are in the minimal surface patch regions, which have typical radius of order 
L/4. As the genus increases, the number of cylindrical handles increases, as will the difference between (i?) and L/4. 
This, in turn, will increase the overall compression energy. 

For our proposed smectic structure to be stable against the uniform flat phase, the positive energy contributions 
from Fc, Fb, and -Fcore must be compensated by a large negative saddle-splay energy. Since there is no chirality in 
our construction the unit cell length can only arise as a result of the different scalings of these energies. We found 
that the scaling competition between the core and saddle-splay energies was capable of generating a length scale of 
L = 50ao, giving a preferred cell size for the P and I-Wp surfaces on the order of: 

^ ^ L 127r|g24| + [16A-6]£ 
L ^ 72|j^24|-2.8£ 

"^^-^p - ^ - im^A-e ■ 

We compare the total energy of the minimal surface smectics against that of uniform smectic configuration. Fig. 0] 
shows the exact stability diagram of both the P and I-Wp surfaces, assuming a crystal unit cell length L ~ 50ao to 
make contact with experiment. There are stable solutions of the P surface smectic for values of \K24\/Ki > 6.5 and 
for the I-Wp surface smectic for values of \K24\/Ki > 2.7. 



B. Numerical Minimization 



In order to study elastic relaxation away from our simplified construction, we performed a numerical minimization 
of the smectic energetics of Eq. and Q on a 3-dimensional grid. The smectic field was represented by a single value 
of the phase <i> at each point on a 41 or 128 unit cubic grid with periodic boundary conditions. Discretized energy 
expressions were used to calculate the local values of compression, curvature, and saddle-splay energy. We manually 
put line defects into our lattice which forced the P surface or I-Wp surface topology, similar to the analysis in 1^. 
The energy was then minimized by the conjugate gradient method. Because of the troublesome nature of defining 
and allowing defects in such a phase field simulation, as well as the relatively small grid size, the numerical values in 
these simulations cannot be considered precise. However, study of the numerical minimal energy configurations are 
instructive for both validating and expanding upon our models. 

We imposed parallel boundary conditions on the smectic field at the defect lines (in other words, the defect lines 
become surfaces of constant phase) . This corresponds to our "small core" construction. The results are pictured in 
Fig- El The shading indicates local bending energy density. For B/Ki relatively small, the middle surface between 
inner and outer defect lattices is close to a minimal surface. For layers away from the middle surface, such as that 
pictured in Fig.^Ia), the curvature energy is concentrated around the cylindrical regions, while the remaining areas 
are close to minimal surfaces. Fig.EJb) shows the strong Fourier components in the — Q plane for the TWp surface 




FIG. 5: Numerical minimization results with parallel defect boundary conditions. Images (a) and (c) show two different smectic 
layers for a relaxed configuration of the I-Wp surface. Surface coloring is proportional to local curvature energy density, with 
blue and violet denoting the highest values. Plot (b) shows the strongest Fourier components in the kz = plane for this 
smectic configuration. 



smectic. The numerically relaxed configurations for moderate values of A = a/ Ki/B m 1 and i/og ~ 40 layers per 
unit cell showed the same off-peak increase for \kmax{0^ 4')\ as was found experimentally by Pansu et al. |6|. 

For higher values of B/Ki, such as in Fig.jSJc), the smoothly curved regions of our construction start to become 
faceted, with the curvature condensing onto discrete folds. This is not unlike the condensation of curvature found 
in crumpled elastic sheets p^ . Motivated by this apparent tradeoff between curvature and compression energies, we 
examined the model for relaxation pictured in Fig.0 in which the regions of high curvature become more concentrated 
in space while the remaining overall surface becomes more polyhedral. This model has one continuous parameter 
X G [0, 1] which interpolates between the minimal surface model at x = 1 and a faceted surface with uniform layer 
spacing away from tilt grain boundaries at x = 0. In this model, we shrink patches of the minimal surface onto vertices 
of a polyhedron with the same symmetry. The patches are connected by cylindrical sections, and the remaining area 
is filled with flat facets. An expression for the compression energy of this structure, using our earlier techniques, is 
complicated, but numerical evaluation for the P surface found it was well approximated by a simple linear function 
of X- The curvature energy is found by adding up the contributions of the cylindrical regions. We have ignored the 
complicated non-linear effects at low x, so our analysis is valid near x — ^- The resulting combined compression and 
curvature energy density for the P surface smectic is 



Fb+Fc ^ 4.32 X IQ-^xBL^ + GttKiL 



Minimizing in x gives 



-log(L/4pe) + (4pjL-l) 



(19) 



X 



4.32 X 



0.94- 



(20) 



where the last equality is for L — bOao. Since A > ao and this treatment is only valid for x ^ Ij we see that curvature 
condensation will not happen for our small smectic repeat unit, but might occur for L ^ 50ao. 



IV. LARGE CORE MODEL: CALCULATION OF CURVATURE ENERGIES USING WEIERSTRASS 

REPRESENTATION 



A. Geometric Construction 



The large core variant of our construction begins with a single continuous surface, then adds successive layers which 
are equally spaced along the original surface normals. As we showed in Section^ the continuation of uniform layer 
spacing will eventually lead to singularities. Fig.[7|shows the results of this construction for the interior of a P surface 
section. The initial curvature singularities occurs on the nearly circular interface between unit cells of the P surface: 
essentially, the circular interface shrinks to a point, pinching off the connection between the interior layers and the 
outside of the unit cell. Any further layers within the pinched layer must be topologically like a sphere, with genus 0. 
Note that this layering scheme appears similar to the construction of constant mean curvature (CMC) surfaces 






FIG. 7; Construction of smectic blue phase based on uniform layer spacing of the Schwartz P surface. The image shows a 
cutaway of the layered structure inside the P surface. 



but is essentially different; there is no energetic preference for constant mean curvature and the CMC surfaces do not 
minimize the compression energy. We shall label the normal distance from the initial minimal surface to the last fully 
connected surface as amax- Since the principal radii of curvature are equal and opposite for points on the minimal 
surface, the value of |amaa;| will be the same on both sides of the minimal surface. We call the region within the 
pinch-off surface, at distance greater than amax from the minimal surface, the core region since our initial ordering 
is disrupted in this volume. Outside the core, the saddle-splay energy per unit cell is Fss — 16|-?i'24|(5 — ^)0'max- 
The values of g and amax for the P and TWp surfaces are given in Table Hll Fig.|3a), (b) and (c) show the surfaces 
defined by a^ax for both minimal surfaces. Within the pictured surfaces, the continuation of our initial construction 
produces additional curvature singularities and lines of surface self-intersection. At the end of this section we justify 
a new model for the core, which fills it with even layers joined by tilt grain boundaries. 

By insisting on uniform spacing away from the fiducial minimal surface, we determine the exact shape of all the 
layers outside the core with no free parameters. Clearly, the shape equations of all the layers can be obtained from 
that of the original surface if we know how all quantities evolve with normal displacement. The curvature evolution 
was derived in Section ^ We can also map the area element of the initial surface to that on successive surfaces by 
dA' = dA X (1 -I- 2aH + a^Kc)- This allows us to express the curvature component of the smectic free energy (PJ as 
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TABLE II: Surface geometries and energies in the "large core" construction for the P and I-Wp surfaces. The first column 
gives the generating function for the Weierstraf5 representation, followed by the maximum half-thickness of continuous layering, 
the curvature energy, saddle-splay energy, and grain boundary core energy for this ansatz. 

where the area integral is taken over the center surface. In our case the area integral is over the minimal surface repeat 
cell and H = 0. We evaluated the surface integral in Eq. H21|) analytically using the Weierstrafi representation for the 
minimal surfaces |l9l l20j , as described in Appendix ^ Table ^ gives numerical values for the bending energy using 
Eq. H21|l and the saddle-splay energies using the Gauss-Bonnet theorem and Eq. ||2J). In both cases the integration 
limit a^nax corresponds to the value of a for which H' in Eq. |(7J) first diverges at some point on the surface. Thus 

ttmax = {—Kq"^) with Kq^" the minimum (most negative) value of Kg on the minimal surface. 

All that remains is to find the optimal way to fill the " core" volumes for \a\ > \amax\- The shape of these regions 
is shown in Fig.[2Ia), (b), and (c). In the core, the continuation of uniform layer spacing from the minimal surface 
leads to curvature singularities and layer self intersections. However, as shown in Fig. [2Ia), we can fill this volume 
with uniformly spaced domains of smectic order which intersect each other in tilt grain boundaries. Within a domain, 
each point on a given layer is an equal normal distance away from some point on the initial minimal surface. Points 
on the grain boundaries between domains are an equal normal distance away from at least two different points on the 
minimal surface, so the grain boundaries occur where the uniform spacing used outside the core would cause layer 
self-intersection. We have no general principle of why the singularities in the core region are domain walls instead of, 
say, focal conic domains (such as you may find by developing a catenoid) . For both the P and I-Wp surfaces uniform 
spacing leads to disclination walls, as it probably does for all but highly symmetric surfaces. 

Thus, with the aim of avoiding layer compression, we choose this solution for the core structure, which should 
serve as an energetic upper bound. The energy of the core will arise from the regular curvature and saddle-splay 
energies within the smectic domains, along with a surface energy at the grain boundaries and a separate line energy 
in the highly distorted regions where the grain boundaries meet. A conservation estimate of the energy of a tilt grain 
boundary is approximately KiUq^ per unit area [2l| (assuming a melted wall of thickness oq at the defect plane). 
Numerically, the total area of grain boundary walls per unit cell in this construction is ^ = 1.65L^ for the P surface 
and A ~ 1.36L^ for the I-Wp surface; the corresponding energies are given in Table ITU The lines at the intersection of 
grain boundaries are located where the -1-1 and +1/2 disclinations were in the small core model of Section llTTl Around 
these lines the director field varies so rapidly that the smectic order should melt completely, leaving a core of radius 
^ flo and line tension e. The mean curvature energy in the smectic domains is relatively small, since the shape of the 
innermost layers is mapped by Eq. Q from the region around the umbilics (flat points) on the minimal surface. The 
inner layers are thus nearly flat. We found that the curvature energy was negligible by numerical evaluation. 

The saddle-splay energy in the core merits special attention. Since the smectic layers in the core are not closed 
surfaces, our simple topological arguments do not apply. We could insist on connecting the surface sections together 
in a natural way across grain boundaries to obtain closed surfaces with the topology of a sphere. However, these 
closed surfaces are nearly polygonal, with all the Gaussian curvature concentrated at the vertices, or along the lines in 
three dimensions where the smectic order has melted. Thus, it is natural to assume that the smectic order parameter 
vanishes in these line-like cores, and the total saddle-splay of the core is nearly zero. 

The dominant free energy terms which determine the stability and preferred size of this construction are therefore 
the bend, saddle-splay and grain-boundary energies. Minimizing these energies in a unit cell of length L gives: 

L = {\Fss\-FB)/i2FGB)- (22) 

At this length the phase is stable against the standard, flat smectic when \Fss\ > Fb- Reading values from Table ITU 
the minimum for the P surface structure occurs for L — 7.06ao(|-ftr24| — 0.523Ki) / Ki, while that for the I-Wp occurs 
at L = 16.6ao(| A'241 — 0.325Ki)/ Ki. To compare with experiment we set L — 50ao, where ao is the smectic layer 
spacing. This requires \K24\/Ki = 6.6 for the P surface or 2.6 for the I-Wp surface. Thus the I-Wp ansatz should be 
the stable phase at any value of if 24- 
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FIG. 8: Numerical minimization results with arbitrary defect boundary conditions. 
P-surface smectic. 



The image shows the core region of a 



(a), 



1(b) r 



^ ■ 



FIG. 9: Discrete Fourier transform of proposed smectic structures. Plot (a) is for the P surfaced based construction, (b) and 
(c) are for the I-Wp construction. Plots (a) and (b) show the largest Fourier components in the positive coordinate octant of 
k space. Plot (c) shows the largest Fourier components in the positive coordinate quadrant of the kz — plane. Point size is 
proportional to the Fourier coefficient at that k. 



B. Numerical Minimization 



As before, we used conjugate-gradient minimization to explore the possible smectic configurations. In this case we 
relaxed the parallel boundary conditions, letting the interface of smectic field and defect line be arbitrary. The results 
are shown in Fig.|Sl The core region resembles our "large core" model, except that the defect walls have disappeared 
at the expense of greater overall compression. Locally, the core layers look like focal conies around the defect lines. 
This result undoubtedly points to configurations close to our model but with lower energy, and most likely stable for 
even lower values of |i^24|/^i- Due to numerical difficulties, these phases could only be made stable for large values 
of L/ag > 80, so we could not compare configurations in the region of interest. 



V. COMPARISON TO SPECTROSCOPIC DATA 



Pansu et.al. 0,0,0 found through X-ray diffraction that the cubic BPsmA^ phase showed the strongest diffraction 
along three mutually orthogonal axes. They also found that the peak wavelength of diffraction decreased away from 
these special directions. In order to compare our model to their data, we have calculated the expected scattering from 
a smectic blue phase built by our "large core" construction. Fig.|51shows the fc-space location of the strongest Fourier 
peaks for both a P surface smectic and an I-Wp surface smectic. In each case, we sampled the smectic phase on a 
cubic grid with 128 points per side. We took the smectic phase to be 27r x 37 x Xn/L, where Xn is the distance of the 
sample point normal to the minimal surface. We chose this number of layers to avoid any commensuration effects. 
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FIG. 10: One patch being uniformly translated along its normal. The top and bottom surfaces along with the perpendicular 
surfaces generated by the normals at the boundaries makes an integration volume V. 

Because there are many layers, this was numerically challenging and a higher precision study is needed, such as the 
study of scattering from triply-periodic minimal surfaces by Garstecki and Holyst . 

For both structures, and we expect for general minimal surface smectics, the Fourier peaks were strongest along 
the umbilics of the minimal surfaces, i.e. normal to the points at which the surfaces are flat. The Fourier peaks of 
the I-Wp structure are a clear match to the observed X-ray diffraction of BPsmA^- However, as Fig. Efc) shows, 
there is no clear deviation in diffraction maximum wavelength as a function of angle (Pansu et.al. reported a 7% 
deviation over 25° angle). This indicates the need for further refinement to our model if it is to explain the current 
data exactly. The discrepancy in off-peak behavior between our simulations and real data could arise from some 
relaxation of our proposed structure to better accommodate the competing energies. As we noted in Section flll Bl the 
numerical relaxation of the I- Wp smectic with parallel ( "small core" ) boundary conditions does produce the correct 
off-peak behavior for the diffraction maxima. 

VI. CONCLUSION 

We have presented a ne w mod el for smectic blue phases which is stable for ph ysic ally realizable values of K24 and 
penetration length A = y^Ki/B. This work refines that of our previous paper im and gives further weight to our 
proposed organizing principle of smectics built on minimal surfaces. Unlike the traditional blue phases, our model does 
not rely on molecular chirality. It would be interesting to add chirality into this model to see how the new length scale 
alters the equilibrium structures. We have also addressed the broader question of geometrical frustration between 
Gaussian curvature and smectic order. Our computed structure factor matches the scattering from the Suibp^ phase 
and so we remain optimistic that our construct could be verified, perhaps through freeze-fracture. 
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APPENDIX A: ALTERNATE DERIVATION OF CURVATURE EVOLUTION: LAYERING AND 
CURVATURE FRUSTRATION IN CURVED SPACE 

The first result in Eq. (|SJ also follows by considering the derivative of the mean curvature along the layer normal: 

= iv.[N(V.N)]-i(V.Nf (Al) 

When the layers are built by developing surfaces parallel to a fiducial surface, the layer normals do not change - in 
other words, since we are translating the surface parallel to N, we must have (N • V) N = 23]. Adding this term 
to Eq. IjAljl we have: 

^ = -V.[N(V.N)-(N.V)N]--(V.N)^ 

= Kg-2H^ (A2) 



dH 
da 



(N-V) 



1, 



-V -N 
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Similarly, consider a patch on one layer and the corresponding patch on a uniformly translated surface as shown in 
Fig- El Because the principle directions of the two surfaces are unchanged, the geodesic curvature of the corresponding 
boundaries are identical. It follows from the Gauss-Bonnet theorem that /^^^ ^ KcdA is the same on the two patches. 

Thus, over the volume V swept out by making the translated patch (shown in Fig. I10|l . Jy ■ dA = 0. Therefore 

V • (i^cN) = and 

= N ■VKg + KgV -N 

= ^ + 2KgH (A3) 
oa 

The intrinsic frustration between curvature and uniform layer spacing can be relieved by considering layered systems 
in curved space, just as double-twist in the classical blue phases can fill the surface of the three-dimensional sphere 
without defects [2J|. It is straightforward to see this using the derivations above. In curved space, we simply replace 
derivatives with covariant derivatives, so that V ^ D where D is the covariant derivative |2^. The mean curvature 
is -ff = ^DiN^ and is essentially unchanged. However, the relation between Gaussian curvature and saddle-splay is 
more subtle. Recall that in flat space, the coefficients of the two terms in the saddle-splay are constrained so that the 
saddle-splay only depends on first derivatives of N. We shall see that this is spoiled in curved space. We have: 

D • [(N • D) N - N (D • N)] = A [NW.jN' ~ N'DjN^] 

= {D.NWjN' - D.N^DjN^) + [Di.Dj] Ni (A4) 



The first term in (|A4|) is ~2Kg, a factor times the Gaussian curvature. Note that in flat space, the covariant 
derivatives become simple derivatives and commute. In curved space this is no longer true and the commutator term 
is just Rnn-, the component of the Ricci tensor in the normal-normal direction. The evolution equation for Kg is 
unchanged since Stoke's theorem (properly modified) holds in curved space. Thus in curved space we have 



dKc 
da 
dH 



-2KgH 



„ = Kg- 2H' - -Rnn (A5) 
oa I 

Thus we can have uniformly spaced layers with H — and Kg constant {Kg may vary in the layer, but will not 
change from layer to layer) if we embed the smectic in a space with Rnn — "^Kg- Since we are considering surfaces 
with Kg < 0, this suggests a space with negative curvature. 

To see this, we consider a special coordinate system in a 3-dimensional layered structure. For unbroken layers we 
can uniquely define the continuous variable $ which labels the layers. Furthermore, we can use $ as a local coordinate 
to define a coordinate basis in which one basis vector 63 is dual to the directional derivative in $, and the other basis 
vectors ei and 62 lie in planes of constant This coordinate system is known as Gaussian normal coordinates 25j. 
Furthermore, for uniform layer spacing the measurement of distance along 63 cannot depend on coordinates 1 and 2. 
This, plus the orthogonality of 63 lets us scale our coordinates such that g^i = 63 • ei = Sat. 

Since we are using a coordinate basis, the covariant derivative is defined with the typical connection coefficients 
r*jj, . Elementary considerations relate the connection coefficients to the extrinsic curvature tensors of the layers taken 
as 2-sheets: 

Ktj = Tsij = -Ti3j = -Tj3i for i,j ^3 

T^33 = r3.3 = r33, = 0. (A6) 

Laborious but autonomic calculations show that the Ricci scalar R — 2Rnn and so the space has negative scalar 
curvature. This is not surprising: in the classical blue phases the saddle-splay elastic constant needed to be positive, 
favoring positive Gaussian curvature, and the resulting structure could be defect-free in positively curved space. We 
have merely "flipped the signs" on the last sentence. Further work on flattening the hyperbolic space into R"^ would 
be interesting and may shed light on the preferred lattices. 



APPENDIX B: WEIERSTASS REPRESENTATION OF MINIMAL SURFACE 



The Weierstrafi representation of minimal surfaces relies on the fact that there is locally a one to one mapping from 
a minimal surface onto the unit sphere via the surface normal pj| . The unit sphere can, in turn, be mapped onto 
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FIG. 11: Weierstrafi representation of the P surface, (a), and I-Wp surface, (b). The shaded region in each is 1/48 part of the 
pictured surface repeat unit and contains all the symmetry of the entire surface. The points represent the location of a; = for 
the Weierstrafi mappings given in Table ITTl 



the complex plane by stereographic projection, with the plane passing through the equator of the sphere. The latter 
mapping gives the components of the surface normal in terms of the complex variable to = a + ir as 



25R uj 



No 



2Q LU 



OJl 



1 



1 



(Bl) 



The reverse mapping between the complex plane and the minimal surface itself can be expressed via a single generating 
function II(lu,lu). In terms of R(a;, w), displacements on the minimal surface are 



dx^ = [(1 - u^) R(cj, Ld)duj] 
dx^ ^ -3 [(1 + uj^) R{lu, Ld)duj] 
dx^ = -2^[ujR{uj,u})dLj] 



(B2) 



The relative coordinates of points on the surface can be found by integrating Eq. (jB2p . 

To complete our reparameterization of Eq. (|21() in terms of we quote expressions for other quantities of interest 
on the surface [l9ll20| : 



l^YdadT 



dA=\R{uj,oj)f (l- 



(B3) 



The Weierstrafi representations for the surfaces we consider are given in Tabled] For the P surface, a; = (z = — 1 
on the unit sphere) corresponds to the point indicated in Fig. Illf a). The entire shaded region in Fig. Illf a). which 
is 1/48 of the P surface repeat cell, maps to the region on the unit sphere bounded by the intersection of the sphere 
with the planes x = 0, y — 0, z = —x and z — —y. For the TWp surface, the position of w = is shown in Fig. llir b'). 
The shaded region is again 1/48 of the 1-Wp repeat cell and corresponds to the region on the unit sphere bounded 
by the planes x = 0, y = 0, z = x, and z = \/2y — x. 
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